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Abstract. It is speculated that the accretion of material onto young protostars is 
episodic. We present a computational method to include the effects of episodic ac- 
cretion in radiation hydrodynamic simulations of star formation. We find that during 
accretion events protostars are "switched on", heating and stabilising the discs around 
them. However, these events typically last only a few hundred years, whereas the in- 
tervals in between them may last for a few thousand years. During these intervals the 
protostars are effectively "switched off', allowing gravitational instabilities to develop 
in their discs and induce fragmentation. Thus, episodic accretion promotes disc frag- 
mentation, enabling the formation of low-mass stars, brown dwarfs and planetary-mass 
objects. The frequency and the duration of episodic accretion events may be responsible 
for the low-mass end of the IMF, i.e. for more than 60% of all stars. 



1. Introduction 

Young protostars grow in mass mainly by accre ting material from their su r rounding 



discs . The accretion is speculated to be episodic (IHerbig]|1 977: Dopita 19781 : iReipurth 



19891 : iKenvon et alJll990fc lEvans et al.ll2009l : iDunham et alJl2010h7 as e.g. in FU Ori- 
type outbursts in which a large fraction o f the star's mass is delivered within a few 
hundred years (Hartmann & Kenvon| [l996l : IZhu et al. 1 12010b . The estimated accretion 



rates during these bursts are up to ~ 5 x 10 M©yr , so each burst can deliver a 
few 10~ 2 Mq onto the protostar. Further evidence fo r episodic accre tion comes from 
the periodically spaced knots seen in bipolar jets (e.g. Reipurth ll989h . Bipolar jets are 
driven off along the protostellar rotation axis by the energy released as material spirals 
in and accretes onto the protostar; periodically spaced knots therefore imply episodic 
accretion. Finally, the luminosity problem provides indirect observational support for 
episodic accretion. By the end of the Class phase (~ 10 5 yr), half the protostar's final 
mass has been accumulated. Thus for a final mass of 1 M Q , the mean accretion rate onto 
the protostar must be ~ 5 x 10^ 6 Mq yr , and the mean accretion luminosity must be 
~ 25 Lq. This is much larger than th e observed bolomet ric luminosities of typical solar- 
type protostars ( Kenyon et al.| [l990: Ev ans et alj|2009h . Episodic accretion mitigates 
this problem, as the lum inosity is only large intermittently, and mo st protostars are 
observed between bursts (Dunh am et all20ld : IOfm er & McKee1 l20Tlh . 
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Accretion of material onto a young protostar results in large amounts of accre- 
tion luminosity which may effect the evolution of the protostar's surrounding disc 
dKrumholzll2006l:lBatel2009i : lQffner et aDl2009l : IUrban et al.ll2010l : lKrumholz et all2010l : 
Offne r et al. 2010h . This is particularly important for the low-mass end of the IMF (e.g. 



for low-mass stars and brown dwarfs) as gravitational fragmentation of massive ex- 



tended discs is considered to be one of the main ways of p roducin 
Whitworth & Stamatellosl 12006b IStamatellos & Whitworthll2009alf 



these objects (e. 



Bolev et al. 201 



Kratter & Murrav-Clay 2011). Thus, in order to determine the full effect of radiative 



feedback in low-mass star formation, it is important to describe in detail the role of 
episodic accretion. 



2. An episodic accretion model for hydrodynamic simulations 

Material spirals from the outer disc region towards the young protostar due to the work 
of the gravitational instabilities in transporting angular momentum outwards in the disc. 
Once this material reaches within a few AU from the protostar the GIs becomes inef- 
fective as the disc is hotter there. The material is then deposited on a notional inner 
accretion disc (IAD), where it piles up until it becomes hot enough that thermal ion- 
isation couples the matter to the magnetic field. At this point the magneto-rotational 
instability (MRI) is activated, transporting angular momentum outwards, and thereby 
allowing the matter accumulated in the IAD to spiral inwards and onto the central pro- 
tostar. 

Simulations of self-gravitating hydrodynamics on the scale of molecular cloud 
cores (i.e. sizes from 10 4 to 10 5 AU and masses from 0.3 to 10 M Q ) can achieve suf- 
ficient resolution to capture the formation of discs around young protostars, and the 
effects of GIs in the outer regions of such discs. However, such simulations cannot 
capture what happens in the inner disc reg ion, where sinks are invoked. Contrary to 
previous studies (e.g. lAttwood et al.ll200"9f) . we assume that the matter which enters a 



sink is not instantly accreted onto to prot ostar but it is acc umulated on the IAD, within 
the sink radius. We adapt the model of IZhu et~ail ( 2010h to phenomenologically de- 



scribe the properties of episodic accretion. The details of the model are discussed in 



Stamatellos et al] (1201 lbl) . An episodic accretion event is initiated when 

Mad > Mmri (1) 

which is presumed to correspond to an IAD temperature of Tm = 1400 K. In the above 
equation 



M^O.lSMel-^- LJ*fr _, 1 . (2) 



f I M lAD ^ 
,0.2M o j ylO^Moyr- 1 

is the mass accreted onto the protostar in an outburst event. The duration of the event 
is set as 

At =* 250 vr (—f 1 ( Mmri \ (3) 
mri Z3U y r ^ .lJ v 0.13M G j' {> 

where q-mri is the viscosity parameter provided by the MRI. To match observations of 
FU Ori stars, we assume an exponential decay of the accretion rate 

M M ri 

M*, EA = — e a 'mri , to < t < t + At MRl . (4) 

AfMRI 
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Once the accretion rate of the protostar has been determined using the above 
model, its luminosity at any time is given by 

/ M \ 3 fGMM 

where M is the mass of the protostar, R its radius, and M is the accretion rate onto 
it. The first term on the righthand side is the intrinsic luminosity of the protostar (due 
to contraction and nuclear reactions in its interior). The second term is the accretion 
luminosity. / = 0.75 is the fraction of the accretion energy that is radiated away at 
the photosphere of the protostar, rather than being expended driving jets and/or winds 
(Offner et al. Il2009h. We assume R — 3Rq is the typical radius of a young protostar 
dPalla & Stahlej[l993i) . At the initial stages of star formation the accretion luminosity 
dominates over the intrinsic luminosity of the protostar. 

The above model describes macroscopically the effects of episodic accretion, which 
in turn regulates the radiative feedback from the protostar. 



3. The role of episodic accretion in low-mass star formation 

To evaluate the consequences of episodic accretion for disc fragmentation and low- 
mass star formation, we perform radiation hydrodynam ic simulations of a collapsing 



turbulent molecular core. We use the SPH code SEREN (Hubbe r et al.ll201 lb. in which 



the ra diative transfer is treated with the diffusion approximation of IStamatellos et al. 



( 2007bh . The outer envelope of the core extends to 5 x 10 4 AU, and the total core 



mass is 5.4 M . We perform three simulations, all with the same initi al conditions, and 



differing only in their treatments of the luminosities of protostars (see IStamatellos et al. 



201 lb, for details). In all three simulations the initial collapse leads to the formation of 
a primary protostar (i.e. a first sink) at t ~ 77 kyr, and this quickly acquires an extended 
accretion disc (Figs.[T|[3]>. The simulations only diverge after this juncture. 



No radiative feedback from protostars. In the first simulation there is no radiative 
feedback from the protostar, and the gas in the disc is only heated by compression and 
by viscous dissipation in shocks. The disc stays cool and fragments resulting in the 
formation of 7 secondary protostars, with masses from 0.008 M Q to 0.24 M . This is 
illustrated in Fig. [T] 



Continuous accretion/continuous radiative feedback from protostars. In the sec- 
ond simulation, we assume that the matter entering a sink is immediately accreted onto 
the protostar at its centre. The accretion of material onto the protostar is continuous 
(see. Fig. HJ top, blue line). This results in an accretion luminosity which is typically 
10 to 100 L Q (see Fig. 01 bottom, blue line). It therefore heats the surrounding disc, 
entirely suppressing disc fragmentation, at all radii, as it is illustrated in Fig.|2] 



Episodic accretion/episodic radiative feedback from protostars. In the third sim- 
ulation, the accretion rate into the protostar during the MRI-driven accretion bursts is 
very high (see Fig. |4j top, black line), resulting in high luminosity (see Fig. |U bottom, 
black line). This heats and stabilises the disc against fragmentation. However, as the 
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Figure 1 . Evolution of the accretion disc around the primary protostar forming 
in a collapsing turbulent molecular cloud core, without radiative feedback from the 
protostar. The disc around the primary protostar increases in mass, becomes grav- 
itationally unstable, and fragments to form 3 low-mass stars, 2 brown dwarfs and 
2 planetary-mass objects. The colour encodes the logarithm of column density, in 
g crrT 2 . 




Figure 2. Evolution of the accretion disc around the primary protostar with con- 
tinuous accretion, continuous radiative feedback from the protostar. The disc grows 
in mass, but radiative feedback makes it so hot that it does not fragment. 




Figure 3. Evolution of the accretion disc around the primary protostar with 
episodic accretion and episodic radiative feedback from the protostar. The disc 
becomes gravitationally unstable (first and second column) but fragmentation is 
damped by heating due to an accretion burst (third column); however, after this burst 
the disc cools sufficiently to undergo gravitational fragmentation; 2 low-mass stars 
form in the disc. 
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Figure 4. The evolution of the accretion rate and the luminosity of primary proto- 
star. Blue lines correspond to the the case with continuous accretion and feedback, 
and black lines for the case with episodic accretion and feedback, (a) shows the ac- 
cretion rate onto the primary protostar amd (b) shows the accretion luminosities, for 
the cases with continuous (blue) and episodic (black) accretion. 



mass of the protostar builds up, the mass that is required to activate MRI also goes up, 
and at the same time the accretion rate into the sink goes down, so that the intervals 
between bursts become longer. During these intervals the luminosity of the protostar 
is relatively low, hence the disc cools down. Within a few kyr of the formation of 
the primary protostar, the interval between successive accretion bursts has increased to 
~ 5 kyr, and this is sufficient time to allow the outer disc at R ~ 50- 150 AU, to undergo 
gravitational fragmentation. This is illustrated on Fig. [3] In the first two frames, at 82 
and 83 kyr, the disc is unstable and tries to fragment, but then at 84 kyr it is heated 
by an accretion outburst and stabilised. Once the outburst is over, the disc cools back 
down fast, and fragments to produce 2 low-mass hydrogen-burning secondaries. 



4. Conclusions 



If episodic accretion is a common pheno menon among young proto s tars as observa- 
tional and theoretical evide nce suggests ( Herbid 19771: iDopital 1 19781: lReipurfhl | 19891 : 
Hartmann & Kenvon| [l99^: iGreene et al.ll2008l : lEvans et al.ll2009l : iPeneva et ailboiol : 
Dunham et alj|2010l) . it may limit the effect of the luminosity of a protostar on its en- 
vironment and create favourable conditions for disc fragmentation to occur. As disc 
fragmentation is predominantly a mechanism for forming low-mass stars and brown- 
dwar fs, episodic accretion may then have a significant influence o n the lower end of the 
IMF ( Sta matellos et al.l2007al : IStamatellos & Whitwortrj|2TO9alJbl : IStamatellos et ali2011al) 
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